Transmission distribution, 'P(lnT), of ID disordered chain: low-T tail. 
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We demonstrate that the tail of transmission distribution through ID dis- 
ordered Anderson chain is a strong function of the correlation radius of the 
random potential, a, even when this radius is much shorter than the de Broglie 
wavelength, kp. The reason is that the correlation radius defines the phase 
volume of the trapping configurations of the random potential, which are re- 
sponsible for the low-T tail. To see this, we perform the averaging over the 
low-T disorder configurations by first introducing a finite lattice spacing ~ a, 
and then demonstrating that the prefactor in the corresponding functional in- 
tegral is exponentially small and depends on a even as a —>■ 0. Moreover, we 
demonstrate that this restriction of the phase volume leads to the dramatic 
change in the shape of the tail of V(lnT) from universal Gaussian in InT to 
a simple exponential (in InT) with exponent depending on a. Severity of the 
phase-volume restriction affects the shape of the low-T disorder configurations 
transforming them from almost periodic (Bragg mirrors) to periodically-sign- 
alternating (loose mirrors). 



I. INTRODUCTION 

All the states in one dimension are localized 
at the scale of a mean free path, l e . This means 
that the typical value of transmission through a 
ID region of a length, £, is T ~ exp(— 2C/l e ). 
Since T is exponentially small, the subject of 
recent theoretical studies 1,2 is the distribution, 
'P(lnT), of the log-transmission (and also viola- 
tion of the "orthodox" ID localization 3 for cer- 
tain correlated disorders 4 ). These studies are 
mainly focused on the body of the distribution 
P(lnT). A separate issue is the question about 
the far tail of the distribution, i.e. the behav- 
ior of P(lnT) at | InT 1 ! > 2C/l e . This ques- 
tion is directly related to a more general con- 
cept of the anomalously localized states in dis- 
ordered conductors 5-8 . In Ref. 9 and in subse- 
quent paper 10 it was asserted that the small-T 
tail is dominated by specific configurations of 
the disorder, V(x), namely, the Bragg mirrors. 
These configurations are illustrated in Fig. 1. 
The potential V(x) = 2V cos(2kpx) opens a gap 
2V centered at energy e = kp/2. The cor- 
responding wave function oscillates with a pe- 



riod n/kp and decays as exp (—jx), where 7 = 
V/(2fcp) -C kp is the decrement. Then we have 
I InT I = 2^C = VC/kp. The important assump- 
tion adopted in Refs. 9, 10 is that, with expo- 
nential accuracy, V (InT) can be found by sub- 
stituting 2V cos (2kpx) into the "white-noise" 
probability, exp — |i e Jjf dxV (x) 2 / Ak F \ , of the 

fluctuation V (x). This yields 9 ' 10 lnP(lnT) = 

^/ e ln 2 T^ /2C. Remarkably, the result coinsides 
with the asymptote of the "exact" solution ob- 
tained by Altshuler and Prigodin 11 using the 
Berezinskii technique 3 . 

The Bragg mirror configurations, V(x) = 
2V cos(2kpx), emerged in Refs. 9, 10 upon apply- 
ing the optimal fluctuation approach 12 ' 13 . This 
approach was specifically designed to deal with 
situations when the result is determined by a 
particular disorder configuration. The above log- 
normal expression for V (T) corresponds to the 
saddle point of the functional integral over dis- 
order configurations. Obviously, the statistical 
weight of an ideal Bragg mirror is zero. Rig- 
orous application of the optimal fluctuation ap- 
proach implies taking into account the configu- 
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rations close to optimal. This procedure corre- 
sponds to the calculation of the prefactor in the 
functional intregral. In most cases 14-17 the pref- 
actor behaves as a power law and, thus, cannot 
compete with the main exponent. 
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FIG. 1. (a) Schematic illustration of the decay of the 
wave function within the Bragg mirror, (b) Solid line: po- 
tential fluctuation corresponding to an "ideal" Bragg mirror; 
Dashed line: "real" Bragg mirror with fluctuating phase, (c) 
Fluctuations of phase result in the fluctuations of position 
of the gap center (dashed line) leaving the width of the gap 
(solid lines) unchanged. 



In the present paper we demonstrate that the 
situation depicted in Fig. la differs drastically 
from Refs. 9, 10 due to a large size of the optimal 
fluctuation. Resulting from this large size, the 
large number of "degrees of freedom" makes the 
prefactor exponentially small, so that, the final 
result for V (In T) is determined by the compe- 



tition of the prefactor and the main exponent. 
More specifically, as illustrated in Fig. lb,c, 
weakly perturbed Bragg mirrors include fluctu- 
ations with phase varying along x. These fluctu- 
ations are "dangerous", in the sense, that they 
result in spatial modulation of the gap center 
(Fig. lc) and, thus, suppress the decrement, 7. 
Large size of a mirror translates into a large sta- 
tistical weight of these dangerous fluctuations, 
i.e. it severely restricts the weight of the effi- 
cient Bragg mirrors. 

As we demonstrate in the present paper, due 
to the reasons listed above, the proper appli- 
cation of the optimal fluctuation approach, i.e. 
taking prefactor into account, has dramatic con- 
sequences for the shape of the tail of P(lnT). 
Namely, 

(i) The log-normal result 9 ' 10 has a "universal" 
form, in the sense, that it contains only the mean 
free path, l t . Thus, it is insensitive to the actual 
value of the correlation radius, a, of the disorder, 
(as long as a <^ l t ). In contrast, we demonstrate 
that, with prefactor taken into account, V (hiT) 
depends on a exponentially strongly even for 
a < l e . 
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FIG. 2. Disorder configuration in which the Bragg mirror 
occupies only a part L op t of the total length C. The decay 
of the envelope of the wave function within the mirror and in 
the rest of the sample is illustrated with dashed and dotted 
lines, respectively. 

(ii) It was assumed in Refs. 9, 10 that the optimal 
Bragg-mirror fluctuation extends over the entire 
region C This is indeed the case for Gaussian 
form of V (InT). However, with V (hiT) having 
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non-Gaussian form, it turns out that the opti- 
mal fluctuation corresponds to the Bragg mirror 
occupying only a part, L opt < C, of the inter- 
val C, as illustrated in Fig. 2. The underlying 
reason for this is that the prefactor makes the 
Bragg mirrors very "costly" . 

The paper is organized as follows. In Sec. II 
we introduce the discretization, which is always 
mandatory for the functional integration. We 
choose the lattice spacing to be finite, which 
is the most convenient discretization for aver- 
aging over disorder configurations of the Bragg- 
mirror type. In Sec. Ill the functional integral, 
which defines P(lnT), is calculated with pref- 
actror in the domain, where the Bragg mirrors 
dominate the low-T configurations. In Sec. IV 
we consider the low-energy domain, where the 
Bragg mirrors, being too costly, become ineffi- 
cient. We demonstrate that relevant low-T dis- 
order configurations in this domain are the loose 
mirrors, which are periodically-sign-alternating 
on-site energies. As we show in Sec. IV, such 
loose mirrors form a well-defined subspace in the 
space of all possible realizations of the on-site 
energies. In particular, they dominate the func- 
tional integral for P(lnT), which we calculate 
with prefactor. In Sec. V we turn to the contin- 
uous limit a — > 0. In contrast to Refs. 9, 10, we 
find that, due to the exponentially small prefac- 
tor in the functional integral, it is loose mirrors, 
extending over a part of the chain, rather than 
the Bragg mirrors, occupying the entire chain 9,10 , 
that dominate V (InT) in this limit. The Sec. VI 
we trace the reason why the loose mirrors are not 
captured in the standard analytical techniques in 
ID. 

II. GENERAL CONSIDERATIONS 
A. Discretization procedure 

To calculate the prefactor of V (InT), it is nec- 
essary, as in any functional integration, to adopt 
some sort of discretization 15 . In this paper we 
simply introduce a finite lattice constant (equal 
to 1), and a finite hopping between the sites 
(equal to 2), so that the problem reduces to 



ID Anderson model. The discrete on-site en- 
ergies, V m , are random numbers; their distribu- 
tion function, P(V m ), has a characteristic scale 
A -C 1, which we identify with r.m.s. 
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The discrete version of the ideal "continuous" 
Bragg mirror V(x) = 2V cos(2/^p x) has the pe- 
riod n and corresponds to the on-site energies 
V m = 2Vcos (2irm/n). Then the discrete analog 
of the energy k F /2 has the form 

t n = 4sin 2 (vr/2n), (2) 

where the energy is measured from the band edge 
(equal to —2). To approximate the continuum, 
unlike Ref. 18, we will focuse on the energy in- 
terval close to the band edge, i.e. n >> 1. On 
the other hand, the energy should be well above 
the fluctuation-tail domain, e < E t , where E t 
is determined from the following consideration. 
As follows from the golden rule, for n 1, the 
mean free path is equal to l e = 8e„/A 2 . Then the 
conductance, G e , can be written as G e = e 1//2 / e . 
The upper boundary of the fluctuation-tail do- 
main is determined by the condition Ge ( ~ 1, 
which yields E t ~ A 4 / 3 . The fact that we con- 
sider energies above E t sets the lower bound for 
the values of n, namely, n <C A~ 2 / 3 . 

Once the discretization procedure is specified, 
the averaging over disorder realizations is well- 
defined. In particular, to calculate the statis- 
tical weight of the Bragg mirrors, providing a 
given value of InT, one has to integrate P(V m ) 
over the deviations of the on-site energies from 
V m = 2Vcos (2nm/n) with a restriction that the 
log-transmission for the set {V m } is fixed and 
equal to InT. Translating the claim made in 
Refs. 9, 10 into the "discrete" language, this 
weight is simply equal to T7 m T(2Vcos [2nm/n]), 
i.e. the deviations of V m from 2V cos (2nm/n), 
that are responsible for the prefactor, can be ne- 
glected within exponential accuracy. Below we 
test this assertion by explicit calculation of the 
prefactor. The result of this test can be summa- 
rized as follows. 

(i) Weakly disturbed Bragg mirrors (see Fig. 3a) 
are indeed the dominating disorder configura- 
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tions, providing a given value of InT, only above 



certain energy, E B ; 
illustrated in Fig. 4. 
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(ii) Even for energies bigger than the pref- 
actor is exponentially small. Whether or not it 
competes with the main exponent 9 ' 10 depends on 
the length, £, of disordered region. 
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FIG. 3. (a) Weakly disturbed Bragg mirror on a lattice; 
8V is the tolerance in the on-site energies, (b) Loose mirror 
with "rigidity" in signs of the on-site energies. 
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FIG. 4. Left: schematic plot of the ID density of states 
smeared by disorder. For energies e < Et the states are 
strongly localized. Right: low-T disorder configurations have 
the form Fig. 3b within the energy domain E t < e < Eb and 
the form Fig. 3a for energies e > Eb- 



(iii) Our most important finding is that within a 



parametrically wide energy domain, Eb > e n > 
E t , the low-T disorder configurations are domi- 
nated by the novel entity, which we call "loose 
mirrors". They are illustrated in Fig. 3b and 
represent the alternating regions of equal length 
n/2; within each region the values of V m are ei- 
ther random, but strictly positive or random, but 
strictly negative. The "phase volume" of these 
configurations is much bigger than that of the 
Bragg mirrors. On the other hand, for these con- 
figurations, at large n, the sign "rigidity" of V m 
within each half-period is sufficient to provide 
the Bragg reflection. 

B. Optimal length of the Bragg mirror 
for a given length of the chain 

Taking prefactor into account has a dramatic 
effect on the structure of the optimal fluctuation. 
To clarify this point, suppose that the Bragg mir- 
ror extends not over entire distance C, but only 
over the interval L < C, see Fig. 2. Denote with 
T L the transmission of the mirror. Then for the 
transmission of the entire interval £ we have 



lnTl = I lnTcl +2 



C-L 

I 



(3) 



where the second term describes the transmis- 
sion through the region outside the Bragg mir- 
ror (Fig. 2). It is apparent, that the first term 
in Eq. (3) increases with L, whereas the second 
term decreases with L. This suggests the follow- 
ing procedure to determine the optimal length 
of the mirror. Denote with Pi(lnT) the dis- 
tribution function of T^. Then the distribution 
function of the total transmission for a given L 
can be written as V L {InT + 2(£ - L)/l e }. The 
fact that the Bragg mirror has an optimal length 
can be expressed in the form 



\\nV(£,lnT) | = 



mm 

L 



lnV L \ llnTl +2 



C-L 

I 



(4) 



It is seen from Eq. (4) that the calculation of the 
small-T tail of the net transmission of the entire 
interval C reduces to the calculation of the func- 
tion Vl (InT), which is the characteristics of the 
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Bragg mirror only. If now we use for P^(lnT) 
the result 9 ' 10 V L (\nT) = exp [- (l e In 2 t) /2L] , 
which was obtained neglecting the prefactor, 
then the minimum in Eq. (4) would correspond 
to L — C, i.e. to the Bragg mirror extending 
over the entire interval C. Below we demon- 
strate that, once the prefactor is taken into ac- 
count, the true optimal fluctuation corresponds 
to a "short" Bragg mirror, L < C, within a para- 
metrically wide interval of |lnT|. 



III. WEAKLY DISTORTED BRAGG 
MIRRORS 

A. Calculation of the functional integral 

Here we consider the case when the devia- 
tions, 5V m , of the on-site energies, V m , from 
the optimal values, V m = 2V cos (2nm/n), are 
relatively small. For concreteness we choose 
the Gaussian distribution of the on-site ener- 
gies, P(V) = 7r- 1 / 2 A- 1 exp(-\/ 2 /A 2 ). To calcu- 
late the prefactor due to small deviations, 5V m: 
we adopt the assumption that 5V m are homoge- 
neously distributed within a small interval (tol- 
erance) 5V <C A (Fig. 3a). On the one hand, 
this assumption leads to a drastic simplification 
of the calculation. On the other hand, as we will 
see below, it yields an asymptotically correct re- 
sult. 

With homogeneously distributed SV m , the sta- 
tistical weight of distorted Bragg mirror, Vl, can 
be easily expressed through the tolerance 8V 
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In Eq. (5) we have assumed that 5V not only 
smaller than A, but even stronger condition 
5V <C A 2 /V is met. We will check this condition 
below. 

We now incorporate the fluctuations 5V m into 
the log-transmission of the Bragg mirror, InT. 
As it was pointed out above, random shifts, Se^ 
of the gap center reduce the decrement 7 = 



V jlkp = Vn/27r within each period. This is due 
to the local detuning from the Bragg resonance. 
Quantitatively, the reduction of the decrement, 
7, can be expressed as 



7(<Jei) = 7 



1 




(6) 



As a result, instead of 27L in the absence of fluc- 
tuations, the expression for | lnT| modifies to 



lnT| = 2nJ>(<Je i ) 

i 

« 2 7 L -n 7 ]T (^j . 



(7) 



Consider now a given period, i, containing n 
sites. Denote with V$ the on-site energies 
within this period. Then the shift, Sei, of the gap 
center for this period can be expressed through 
via a discrete Fourier transform as follows 



5e, 



ir\2Y n i \/ (i) sin(^ 



n) 



yr 



, \ V cos' 7 "' 



/ nm 
V n 



(8) 



where the summation is performed over the sites 
within the i-th period. Obviously, for an ideal 
Bragg mirror, V m = 2V cos (2nm/n), we obtain 
from Eq. (8) that 5ti = 0. In the presence of fluc- 
tuations, SV m , the typical value of is propor- 
tional to 5V and can be estimated from Eq. (8) 
as follows. The numerator is the sum of n ran- 
dom numbers, each being ~ 5V. Thus, the typ- 
ical value of the numerator is n x l 2 b~V . On the 
other hand, the denominator is equal to nV/2. 
Then we obtain 



Sei 



Ce 2n (5V" 
(2n)V2 ^ V t 



n 2 C f5_V" 

2 l/2 n 5/2 I y 



(9) 



where the constant C is of the order of 1. 

Looking at Eq. (6), it might seem that the 
condition 5V <C V of the weak distortion of the 
Bragg mirror by fluctuations, and the condition 
Sti <C V of the weak reduction of the decrement 
are quite different. It turns out, as we will see 
later, that 5V <^ V insures that \ji — 7] <C 7, 
and thus justifies the expansion of 7(5e«) used in 
Eq. (7). Substituting Eq. (9) into Eq. (7) we get 
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-8tt 3 C 2 



v V 2 
n 4 V 3 



Now it is easy to see that there exists the optimal 
tolerance 



(10) 



SV = 5V opt 



VAn 5/2 n 3/2 A | lnT | 



Using the fact that the r.h.s. in Eq. (10) is much 
smaller than |lnT|, we can express V through 
| In T | as follows 

7i\ In 71 



4vr 2 AttCL ' 

for which | In Vl | is minimal and is equal to 

ln 2 T /7T^ 2 



(13) 



V 



nL 



+ 8ttC 2 



L 3 5V 2 



\^ L \ = ^- L (-) + LA (T) , (14) 



(11) 



n 2 |ln 3 T|' 

Further steps are straightforward. Using Eq. 
(11), we can rewrite the exponent in Eq. (5) 

as 



where 



A (T) = In 



A 



5V, 



In 



4ttCL 



opt , 



( n 3/2| lnT | 



(15) 



1iiPl| = | 
vr 2 ln 2 T 



Lin 



A\ V 2 7 
WJ + 2A 2 



2n 2 A 2 7 



+ Uln 



/ A_\ 8tt 2 C 2 L 3 ( 5U 2 
W/ + n 3 A 2 ln 2 T 



• (12) 



depends on T very weakly. It is also seen from 
Eq. (15) that the C enters into the final result 
only as a factor under the logarithm, so that our 
assumption about homogeneous distribution of 
5V m is justified. 



Now, in order to calculate the tail of the transmission distribution, V (C,\nT), we substitute 
Eq. (14) into Eq. (4) 



|ln ^ l = m i n {2rG0 



In Tl +2 



C-L 

~17~ 



+ LA (T) . 



(16) 



Next we perform minimization with respect to L. This yields the following equation of the optimal 

L = L opt 



7r| In Tl 
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2AnA 



1 + 2 



'C-L 



opt 
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(17) 



Since we are interested in anomalously low transmissions, |lnT| ^> C/l € , the second term in the 
first square bracket in Eq. (17) is small. The second term in the second square bracket contains 
an additional parameter ~ | lnT|/L op t(n 2 A 2 L). Since L = 8e n /A 2 , the combination n 2 A 2 L is ~ 1. 
Thus, the above parameter reduces to \ lnT\/ L opt , which is also small. More precisely, it is of the 
order of l~ x l 2 . Neglecting second terms in both square brackets, and substituting L opt from Eq. (17) 
into Eq. (16), we arrive at the final result 



lnP(lnT,£) 
lnP(lnT, C) 



I 7T\/2A 

l v ^A~ / 

1 / 7T 

2C VnA 
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nAC > | In T| > n 2 A 2 £ 



7T 



) 2 ln 2 T + £A(T) , \\nT\ > (^^j nAC 



(19) 



It is instructive to rewrite the above result in terms of energy, e n ~ (n/n) 2 , and conductance, 
G n » l/(n 3 A 2 ) 
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FIG. 5. Solid line: our main result Eqs. (18), (19) for 
the low-T tail of the transmission distribution. Dashed line: 
log-normal V(T) of Refs. 9, 10. Inset: the portion of the chain 
occupied with Bragg (loose) mirror. 

The behavior Eqs. (18), (19) is illustrated 
schematically in Fig. 5. We see, that the log- 
normal dependence of Refs. 9, 10 takes place only 
for very small transmission coefficients | lnT| > 
C/\fQ. Only in this domain the Bragg mir- 
ror extends over the entire interval, and the 
prefactor [second term in Eq. (19)] is smaller 
than the main exponent. Within a wide do- 
main C/Vk > |lnT| > C/l e the probability, 
V (CAnT), behaves simple exponent. The 
underlying reason for such a behavior is that the 
dependences of the main exponent and of the 
prefactor on L are opposite. As a result, there 
exists an optimal mirror length 



J opt 



7r|lnT| 



2AnA 



(22) 



which leads to the result Eq. (18). Simple ex- 
ponent Eq. (18) corresponds to the situation 
when L opt is shorter than the length of the in- 
terval, C. More precisely, the portion of the in- 
terval C, occupied by the Bragg mirror, is given 
by Lopt/C = 7T |mT| /V^AnAC. Within the do- 
main C/\^T e > | lnT| > C/l e this portion changes 



from 1 to a small value (A/G) 1 ^ 3 < 1, as illus- 
trated in the inset in Fig. 5. 

B. Justifications of the assumptions 

The above calculation was based on three as- 
sumptions 

(i) 5V <C A 2 /V; we used this condition in the 
expression Eq. (5) for the probability, Vl- 

(ii) 5V <C V; this is the condition that the Bragg 
mirror is well defined. It was also used in deriv- 
ing Eq. (5). 

(iii) Sei <C V; this condition was used in expan- 
sion Eq. (7). 

From the result Eq. (13) of the above calcula- 
tion we find SVopt/V = An 5 / 2 /(47r 2 ). Thus, the 
assumption (ii) is valid under the condition 

An 5 / 2 < 1. (23) 

Below we demonstrate that the same condition 
Eq. (23) guarantees the validity of the other two 
assumptions. 

Assumption (i). Within the domain C/^/T e > 
|lnT| > C/l e , where V (C, InT) behaves as a 
simple exponent, the amplitude of the Bragg 
mirror, given by V = n\\nT\/nL opt , is equal 
to V = v / 2AA and does not depend on £, 
where, as follows from Eq. (15), A = A (L opt ) = 

In (47r 2 C/v / 2AAn 5 / 2 ) « In (l/An 5 / 2 ). To check 

the assumption (i) we rewrite the ratio V5V/A 2 
as (<5VyV)(V 2 /A 2 ) = (An 5 / 2 /47r 2 )(V 2 /A 2 ) = 
An 5 / 2 In (l /An 5 / 2 ) . We see that V < A 2 /5V 

holds under the condition An 5 / 2 <C 1, which is 
precisely the condition (23). 

Assumption (iii) From Eq. (9) we have the fol- 
lowing estimate for the ratio dti/V <C 1 

Se, = (5V\ m 

V 2 1 /2 n 5/2 ^ V 2 J ■ y > 

Substituting into this equation the optimal value 
8V opt = VAn 5 / 2 /(47r 2 ), we obtain 5e/V ~ A/V = 
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[2A(L op t)]- 1 / 2 . On the other hand, A (L opt ) = 
ln^l/An 5 / 2 ^) is large under the same condition 

Eq. (23). This large logarithm justifies the as- 
sumption (iii). 

In conclusion of this Section we would like to 
make the following two remarks: 

(1) The expression for the decrement 7 = 
is the result of the two-wave approximation, 
within which propagating and Bragg-reflected 
waves are coupled only in the first order, i.e. 
by a single harmonics of periodic potential. For 
two-wave approximation to be valid, the second- 
order coupling matrix elements must be much 
smaller than V. The estimate for these second- 
order elements is ~ V 2 /e n ~ V 2 n 2 . Thus, the 
two-wave approximation is valid if V 2 n 2 <C V, 
i.e. Vn 2 < 1. As it is seen from Eq. (23), 
V < A 1/2 n~ 5/4 . Thus Vn 2 < U 1 / 2 ^ 5 / 4 ^ 



An 5/2 



1/2 



n 2 = 



fn 1 ! 2 . We see that the applicabil- 
ity condition of the two-wave approximation is 
weaker than the main condition An 5 / 2 <C 1. 

(2) The condition Eq. (23) implies that the en- 
ergy e n exceeds A 4 / 5 . This, in turn, suggests that 
the conductance G n = kph for e = e n is equal 
to G n = (A 2 n 3 ) , and is large by virtue of this 
condition. 

IV. "LOOSE" MIRRORS 

A. Density of the loose mirrors 

We now turn to the case of low energies. More 
precisely, we consider the domain E B > e n > E t 
(Fig. 4). The upper boundary of this domain 
corresponds to An 5 / 2 ~ 1, whereas the lower 
boundary corresponds to An 3 / 2 ~ 1. For ener- 
gies e n > Eb the transmission is dominated by 
weakly disturbed Bragg mirrors, as discussed in 
the previous Section. For energies e < E t we 
have G e < 1, i.e. these energies correspond to 
the tail states. 

As we enter the low energy (large-n) domain, 
the key component of the above scenario of 
weakly disturbed mirrors gets violated. Namely, 
at n ~ A~ 2 / 5 we have 5V ~ V. This implies 
that almost sinusoidal Bragg mirror cannot re- 
tain its role as an optimal fluctuation, which is 



responsible for low-T values. 

In general, optimal fluctuation constitutes a 
saddle-point in the functional space. In the pre- 
vious Section, by demonstrating that the disor- 
der configurations, contributing to the functional 
integral, differ weakly (by 5V <C V) from the 
optimal configuration, we have justified that the 
saddle point is well defined, or, in other words, 
the expansion around the saddle point yields a 
narrow width of the Gaussian in the functional 
space. In this Section we demonstrate that in 
the energy domain E B > e n > E t there exists 
a well-defined subspace of all realizations of the 
on-site energies, {Kn}, which assumes the role 
of a saddle point. We dub the elements of this 
subspace as "loose" mirrors. A loose mirror is a 
configuration of alternating regions of n random, 
but positive V m and n random, but negative V m . 
It is illustrated in Fig. 3. Obviously, the sta- 
tistical weight of the loose mirrors is small. It is 
easy to see that this weight is equal to 2~ L . Most 
importantly, despite the randomness of V m , the 
fact of the sign rigidity within each interval of 
length, n, is sufficient for the formation of the 
Bragg gap with the well-defined width, and thus, 
for generating the low transmission coefficients. 

The key element of calculation of V(lnT) in 
the regime of weakly distorted Bragg mirrors was 
the expansion Eq. (7), which expressed the fact 
that the decrement 7 weakly fluctuates from pe- 
riod to period. It turns out that in the regime of 
loose mirrors, An 5 / 2 3> 1, these fluctuations are 
also weak. This can be seen from Eq. (24). Since 
in the regime of loose mirrors the only scale for V 
and SV is A, Eq. (24) yields <5q/V ~ (An 5 / 2 )" 1 , 
which is small in the regime of loose mirrors. 
Thus, the expansion Eq. (7) is applicable in this 
regime as well. 

To calculate the distribution V(lnT) in the 
regime of loose mirrors, the calculation in the 
previous Section should be modified in the fol- 
lowing way. For loose mirrors the "period" con- 
sists of interval of n positive V m followed by 
an interval of n negative V m . The magnitude 
of the gap 2V and corresponding decrement, 
7 = nV/2n, are determined by discrete Fourier 
component of this realization of the on-site ener- 
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gies. Then the expansion analogous to Eq. (10) 
takes the form 

|lnT|«— . (25) 



IT 

Then the corresponding expression for V, analo- 
gous to Eq. (11), reads 

v-^Il. (26 > 

It is obvious that the typical value of V is ~ A 
with variance is ~ A/n 1 / 2 . It can be demon- 
strated that the full distribution of V is given 
by 

1/2 

! exp 




(28) 



where the power L/n emerges from the prod- 
uct over periods. With Gaussian p(V), given by 
Eq. (27), this optimization can be performed an- 
alytically in a similar fashion as in the previous 
Section, yielding 



p(V) 



n 



7T&A 2 



n(V-piAf 



&A 2 



7r\ lnTl 



(27) 



J opt 



1/2 > 



(29) 



where the constants Pi and /3 2 depend on the 
actual distribution P(V). For Gaussian P(V) 
they assume the values Pi = Air~ 3 ^ 2 and p 2 = 
1/2 — 1/tt. 

Analogously to Eq. (5), the actual calculations 
reduces to optimization with respect to L of the 
product 



where Pq = p 2 l^.2 + Pf is the constant of the 
order of 1. The corresponding value of the gap 
width is 



vrllnTl 



opt 



nL, 



opt 



= API'\ 



(30) 



The result Eq. (29) is quite similar to Eq. (22), and differs only by replacement of the logarithmic 
factor, 2A, by a constant /3 , which is of the order of 1. Correspondingly, the final results for 
P(lnT,£) are quite similar to Eqs. (20)-(21) 





hxV(\uT.C) . i,,, -| r::..:i::| hr/1 - I ^— ! > lh.i.71 > 

lnP(lnT,£) 

\P 
r 1/2 

where f3 e ff = P — Pi / p 2 - For Gaussian distribution of the on-site energies we have /3 e // 
The results Eqs. (31), (32) were obtained assuming that loose mirrors are well-defined entities, in 
the sense, that the subspace that they constitute within the functional space has a sharp boundary. 
In the next subsection we examine the width of this boundary and demonstrate that this width is 
indeed relatively small. 



(31) 



, (32) 



0.46. 



B. Tolerance of the loose mirrors 

In order to examine to what extent the loose 
mirrors are well defined, we consider below two 
generic sources of violation of the sign rigidity, 
which are illustrated in Fig. 6. 

(i) We allow the on-site energies within "pos- 



itive" periods to assume slightly negative val- 
ues, restricted by — W (see Fig. 6a), and the on- 
site energies within "negative" periods to assume 
small positive values, restricted by W <^ A. 
This allowance increases exponentially the num- 
ber of configurations constituting the loose mir- 
rors. On the other hand, such an allowance sup- 
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presses the gap. As a result of these competing 
trends, there exists an optimal value of W, that 
maximizes V(lnT, £). 

(ii) We allow a small portion, k, of on-site ener- 
gies to assume the "wrong" sign preserving their 
magnitude ~ A. This allowance also increases 
the number of loose mirrors and suppresses the 
gap. Thus there exists an optimal k < 1, which 
we calculate below. 

The quantitative characteristics of the "qual- 
ity" of the loose mirror is the fluctuation, 5e, of 
the gap center due to the above violations, which 
is analogous to the tolerance 5V of Bragg mirror 
in the previous Section. 



are also sign-preserving and ~ A. Thus we have 
5e ~ (n/n)(W/ A). Finally, the product 

2" 




FIG. 6. Possible violations of the sign rigidity of a loose 
mirror: (a) small (~ W <C A) on-site energies with "wrong" 
sign are allowed; (b) sparse "large" (~ A) on-site energies 
having the "wrong" sign (hash-marked lines) are allowed. 

The enhancement of the portion of the loose 
mirrors due to allowance, W, can be estimated 
as (1 + W/A) L ^ ~ exp(WL opt /A). This is an 
exponential "gain" inP(lnT, C). The "loss" due 
to suppression of the gap, similarly to Eq. (12), 
can be expressed as ~ exp [— L opt (Se/ A) 2 }. The 
relation between 5e and W can be established 
from Eq. (8). Indeed, all the terms in numera- 
tor are of the same sign and of the order of W, 
while the corresponding terms in denominator 



cxp 




(33) 



of the gain and loss has a maximum at 
W opt ~ n 2 A 3 . We see, that the allowance 
is relatively small, since W op t/A = n 2 A 2 = 
{l/n)(E t /e n ) 3/2 < 1. This suggest that the 
allowance does not change the result Eq. (31), 
since the correction to ln'P(lnT) due to al- 
lowance, W op t, amounts to a small portion ~ 
W opt / A <C 1 of the main exponent Eq. (31). 

Strictly speaking, the optimal value of the al- 
lowance, W, is well defined if the exponent in 
Eq. (33) is much bigger than one. Upon substi- 
tuting W op t hito Eq. (33), we obtain <5 1 In 7^ | ~ 
n 2 A 2 L op t ~ raA|lnT|. From here we conclude, 
that for large enough |lnT| > l/(nA), the al- 
lowance, W, indeed leads to 5|ln"P| 1. For 
smaller | lnT| < l/(nA) the gain does not play 
a role, so that the allowance, W, is determined 
exclusively by the second term in the exponent 
in Eq. (33). This term falls off at characteristic 
W = nA 2 /L]jp 2 . This allowance is bigger than 
Wopt, since W/W opt w (nA\ lnT])- 1 / 2 > 1, but 
still smaller than A. The latter can be seen if 

1 /2 

we rewrite W/A = nA/L opt in the form W/A 
(W / " opt /A) 3/4 |lnT|- 1 / 2 , which is the product of 
two small numbers. The fact that for small 
|lnT| < l/(nA) the allowance increases with 
| lnT| can be interpreted qualitatively as follows. 
The smaller is | In T"J , the less effort is required 
to create a disorder configuration with low trans- 
mission T. 

The enhancement of the portion of the 
loose mirrors due to allowance, k, is equal 
to ( y 2 K ) Lopt/n = exp(/t(L/n)ln2). The 
loss can be estimated analogously to the 
case (i). Namely, due to the sites with 
wrong sign of on-site energies the ra- 
tio E n m=1 V^sin(^)/E n m=1 V^cos(^) in 
Eq. (8) is ~ k. This yields the estimate 
5e ~ {n/n)K, so that, analogously to (i), the 
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product of gain and loss can be written as 



exp y^ opt 



K 



/ TXK \ ' 



(34) 



This product is maximal for K opt = n A 2 . We 
see that K opt ~ W opt /A, and thus is small, as 
discussed above. 

The fact that the optimal allowances W opt /A 
and K opt are small justifies that loose mirrors are 
well-defined entities. 

V. CONTINUOUS LIMIT 

In this Section we establish the relation be- 
tween the above consideration on the lattice and 
the results of Refs. 9, 10, obtained within the 
continuous approach. To establish this relation, 
we restore the lattice constant, a, in the disper- 
sion law, i.e. e{k)a 2 = 4 sin 2 (ka/2), where k 
is the momentum. For lattice constant a = 1 
the dimensionless parameter that separates the 



regimes of weakly disturbed [Eqs. (18), (19)] 
and loose [Eqs. (31), (32)] mirrors was equal 
to nA 2//5 . To incorporate the arbitrary lat- 
tice constant, it is convenient to first express 
this parameter through the conductance G for 
a — 1. From the relation G n = ??r 3 A~ 2 we find 

A 2/5 = n -3/5 G -l/5_ ThuSj nA 2/5 = ^2/5 j ^1/5. 

For arbitrary a, the number n should be re- 
placed by (ka)' 1 , while G n should be replaced 
by G(k) = {e{k)/E t ) 3/2 , where E t is the up- 
per boundary of the tail states. For a = 1 this 
boundary is expressed through A as E t = A 4//3 . 
Thus, the parameter nA 2 / 5 for arbitrary a trans- 
forms into G{k)- 1 ' b {ka)- 2,b . It is seen that this 
parameter contains the lattice constant and in 
the white noise limit a — * 0, considered in Refs. 
9, 10, it is much bigger than 1, which corresponds 
to the regime of loose mirrors. We thus con- 
clude, that for small a the distribution function 
P(lnT) is given by Eqs. (31,32). With a re- 
stored, these expressions take the form 



lnP(lnT,£) 
\nV{\nT,C) 



0, 



27T 3 G(k) 



ka 



1/2 



7T 3 G 



ln 2 T 



lnTl 



2ir 3 fiG 



1/2 



1/2 



7T 



3/2 



kC 



G^ika) 1 / 2 



> llnTl > 



kC 

nG 



- " (I) f 



lnTl > 



1/2 r 



7T' 



3/2 



kC 



GV^ka) 1 / 2 



(35) 
• (36) 



The result of Refs. 9, 10 correspond to the first 
term of Eq. (36). We see, however, that this 
result, obtained neglecting the prefactor, does 
not survive the white- noise limit a — > 0. For- 
mally, taking the prefactor into account, shifts 
the domain of applicability of the log-normal dis- 
tribution to very small transmission coefficients 
|lnT| > (PokC 2 /-n z Ga) 1 ^ 2 , so that this domain 
vanishes when a — > 0. Physically, the result of 
Refs. 9, 10, does not apply in the white-noise 
limit due to the huge phase fluctuations, that are 
allowed for small a. These fluctuations forbid the 
formation of a Bragg mirror with a weakly dis- 
torted sinusoidal shape of Fig. 3a. Instead, the 
relevant fluctuations have the form of loose mir- 
rors shown in Fig. 3b, where only the positions of 
sign changes are adjusted to the de Broglie wave 



length, 2ir/kF, of the electron. In this regime the 
tail of V (InT) is described by a simple exponent 
Eq. (35) with a non- universal coefficient oc a~ 2 . 
Clearly, in the "continuous language" , the lattice 
constant a should be identified with the small- 
est scale in the problem, namely, the correlation 
length of the random potential. Thus, we arrive 
at the conclusion that the correlation length de- 
termines the coefficient in front of | lnT| in the 
leading term of "P(lnT). In terms of the cor- 
relation length and dimensionless conductance, 
the portion of the sample occupied by the loose 
mirror is given by |lnT| (Ga/ flokC 2 ) 1 ^ 2 . 

Finally, we establish the energy interval, where 
the loose mirrors, and thus Eqs. (35), (36), de- 
termine the far tail of VilnT). For this purpose, 
we equate the parameter G{ky 1 ^ {ka)~ 2 ^ to 1 
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and express the energy Eb in Fig. 4 through the 
correlation length a as E B ~ E t /(E t a 2 ) 2 ^ 5 . This 
yields the sought interval 



E t < e < 



E t 



(E t a 



2^/5' 



(37) 



Recall that E t , the position of the boundary of 
the tail states, does not depend on a. This is 
valid for small a, such that (E t a 2 ) <C 1, i.e. the 
interval (37) is broad. For Eq. (37) to apply 
we should require that at the upper boundary of 
the interval (37) the corresponding momentum, 
k B , is much smaller than the inverse correlation 
length. It is easy to see that this is indeed the 
case, since ksa = aE^ 2 = (E t a 2 ) 3 ^ 10 <C 1. 

Fig. 7 illustrates the main qualitative outcome 
of our consideration. Namely, for a short-range 
potential with a correlation radius, a <C kp 1 , the 
low-T disorder configurations for energies within 
the interval Eq. (37) have the shape of loose mir- 
rors depicted schematically in this figure. 




FIG. 7. "Continuous" realization of a loose mirror for 
short-range disorder with correlation radius a <C kp 1 . 

VI. CONCLUSION 

In conclusion, let us address the relation be- 
tween our results and the analytical results in 
ID, predicting the shape of P(lnT). Neither 
of the "exact" techniques 3 ' 19 allows to pinpoint 
the actual disorder configuration, responsible for 
low-T values. Although they are believed to 
be exact, each of these techniques contain a 
step at which mirror-like configurations are lost. 
Let us illustrate this point using the Berezinskii 
technique 3 as an example. In Fig. 8 a three- 
impurity scattering configuration, employed in 
Ref. 3 (see also 20 ) to make the case for complete 
localization in ID, is depicted. As was explained 



by Berezinskii 3 , the key ingredient of the tech- 
nique Ref. 3 is the observation that the scatter- 
ing paths I and II correspond to the same accu- 
mulated phase 



= 2kp{x2 + %s — 2x\), 



(3* 



and, thus, interfere constructively for any <fi. 
However, within the "exact" technique, the value 
of is assumed to be random, and the averag- 
ing over (f) is performed. Similar procedure is a 
key element of the technique Ref. 19. Calculat- 
ing the higher-order diagrams in Ref. 11 takes 
into account increasingly large number of pos- 
sibilities of constructive interference of different 
paths, all of which are of the type Fig. 8 (corre- 
spond to the same accumulated random phase). 
However, each step involves averaging over this 
phase. In contrast, the Bragg-like configurations 
are those sparse realizations, for which the phase 
accumulated upon traversing the period, 7r/kp, 
first forwards, and then backwards, is not ran- 
dom, but close to 27T. Thus, in our opinion, the 
complete coincidence of the estimate for V(lnT) 
based on the Bragg mirrors and of the result 11 is 
accidental. 

















< 











— 













*3 



FIG. 8. Illustration of the simplest building block of the 
Berezinskii technique Ref. 3. In course of moving along the 
trajectories, I and II, an electron accumulates the same ran- 
dom phase 4> = 2I(f(x2 + x A — 2x\). 

Finally, let us briefly formulate the main mes- 
sage of the present paper. Creating low-T dis- 
order configuration in ID demands this configu- 
ration to possess a long-range order, adjusted 
to the wave vector, kp. Ideal configurations 
with such a long-range order, are of the Bragg- 
mirror type. However, they are very "costly" 
to maintain over a large distance. This is due 
to the phase fluctuations, that tend to violate 
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the Bragg condition. These fluctuations are not 
captured at the stage of calculating the saddle- 
point. They show up at the next stage, i.e. cal- 
culating the prefactor. We have demonstrated 
that loose mirrors, illustrated in Fig. 7, in which 
the long-range order is present, but relaxed, are 
much "cheaper" to create than the Bragg mir- 
rors. On the other hand, as follows from the 
analysis that we have performed, a loose mir- 
ror, shown in Fig. 7, still constitutes an efficient 
low-T disorder configuration. The smaller is the 
correlation radius, a, of the disorder, the wider 
is the energy interval within which loose mirrors 
dominate the low-T tail of the transmission dis- 
tribution. 

Lastly, we are not aware of any numerical work 
in which the tail of P(lnT) was studied close to 
(but well above) the band edge. Recent simula- 
tions are mostly focused at the body of the dis- 
tribution both in ID 1,21 and in 2D 22 ' 23 , and are 
aimed at testing the scaling hypothesis. With 
regard to the tail of the transmission distribu- 
tion, the related characteristics, namely, the den- 
sity of anomalously localized states, was studied 
numerically only in two 24 and three 25-28 dimen- 
sions. 
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